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ABSTRACT 

Existence of planets is binaries with relatively small separations (around 20 AU), such as a Centauri 
or 7 Cephei poses severe challenges to standard planet formation theories. The problem lies in the 
vigorous secular excitation of planetesimal eccentricities at separations of several AU, where some 
of the planets are found, by the massive, eccentric stellar companions. High relative velocities of 
planetesimals preclude their growth in mutual collisions for a wide range of sizes, from below 1 km up 
to several hundred km, resulting in fragmentation barrier to planet formation. Here we show that rapid 
apsidal precession of planetesimal orbits, caused by the gravity of the circumstellar protoplanetary 
disk, acts to strongly reduce eccentricity excitation, lowering planetesimal velocities by an order of 
magnitude or even more at 1 AU. By examining the details of planetesimal dynamics we demonstrate 
that this effect eliminates fragmentation barrier for in-situ growth of planetesimals as small as < 10 
km even at separations as wide as 2.6 AU (semi-major axis of the giant planet in HD 196885), provided 
that the circumstellar protoplanetary disk is relatively massive, ~ 0.1M©. 

Subject headings: planets and satellites: 
binaries: close 



formation — protoplanetary disks — planetary systems 



1. INTRODUCTION. 

About 20% of planets detected via stellar radial ve- 
locity variations reside in binaries (Desidera & Barbieri 
2007). The majority of these systems are wide separa- 
tion binaries, with semi-major axis a b > 30 AU. At the 
same time, four relatively small separation binaries with 
a b « 20 AU (HD 196885, 7 Cephei, Gl 86 and HD 41004; 
Chauvin et al. 2011) are also known to harbor giant plan- 
ets with projected masses M p ismi « (1.6 — i)Mj. In 
these systems the mass of the secondary star (we call 
"secondary" the binary component other than the star 
orbited by the planet, which we denote as "primary") 
M s is found to be close to 0.4M Q and binary eccentric- 
ity e b is close to 0.4. Also, Dumusque et al. (2012) have 
recently announced an Earth-mass companion to a Cen- 
tauri B, a member of the binary (or, possibly, a triple) 
with a b = 17.6 AU, e b = 0.52, and M s = 1.1M . In this 
system planet orbits the star at ~ 0.04 AU separation. 

The uniqueness of these systems lies in the fact that 
forming planets in them is known to provide extreme 
challenge to planet formation theories (Zhou et al. 2012). 
With the exception of a Cen and Gl 86, planets in these 
binaries reside in rather wide orbits, with planetary semi- 
major axes a p i w 1.6- 2.6 AU. In-situ formation of these 
gas giants is expected to proceed through continuous ag- 
glomeration of planetesimals at these locations, starting 
from very small objects (easily < 1 km). However, grav- 
itational perturbations from the eccentric stellar com- 
panion inevitably result in rapid secular evolution (Hep- 
penhcimer 1978), driving planetesimal eccentricities far 
above the level at which bodies can avoid destruction in 
mutual collisions (Thebault et al. 2008). This problem, 
which is often called collisional or fragmentation barrier, 
is especially severe for small planetesimals, 1 — 10 2 km 
in size, for which the ratio of binding to kinetic energy is 
small. It is also more pronounced far from the primary, 
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where the secular forcing by the companion is strongest 
and planetesimal eccentricities are high. 

Marzari & Scholl (2000) suggested that a combination 
of secular forcing by the companion and gas drag act- 
ing on small (1 — 10 km) planetesimals leads to apsidal 
alignment of their orbits, resulting in smaller relative ve- 
locities, and allowing colliding objects to grow. How- 
ever, Thebault et al. (2006, 2008) demonstrated that the 
planetesimal size-dependence of apsidal alignment acts to 
break orbital phasing between objects of different sizes, 
resulting in high velocity collisions between them and 
reinforcing collisional barrier. 

Interestingly, most studies of planetesimal growth in 
small separation binaries have included the effect of the 
protoplanetary disk on planetesimal dynamics only via 
associated gas drag (Thebault et al. 2004, 2006, 2008, 
2009; Paardekooper et al. 2008; Paardekooper & Lein- 
hardt 2010), without accounting for the gravitational ef- 
fect of the disk. Batygin et al. (2011) have considered 
disk gravity in the context of planet formation and evolu- 
tion in systems with highly misaligned, distant (10 2 — 10 3 
AU) stellar companions, affected by the Lidov-Kozai ef- 
fect (Lidov 1961; Kozai 1962). However, this effect is 
probably irrelevant for planetesimal dynamics in small 
separation (tens of AU) binaries, which are likely copla- 
nar with circumstellar disks. 

In this Letter we show that apsidal precession of plan- 
etesimal orbits induced by disk gravity dominates secular 
evolution of planetesimals at separations of several AU. 
As a result, relative velocities at which bodies collide are 
reduced, sometimes by more than an order of magnitude. 
In massive disks this effect presents a natural solution of 
the fragmentation barrier issue for the in-situ formation 
of the giant planets in small separation binaries, such as 
7 Cephei. 

2. SECULAR EVOLUTION. 

We consider planetesimal motion as Keplerian mo- 
tion around the primary perturbed by the gravity of the 
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companion, that moves on larger, eccentric orbit, and 
the disk. Mass of the primary is M p , and we define 
fi = M s /(M p + M s ). We assume eccentricity of the stel- 
lar binary e b to be small and planetesimal orbit to be 
coplanar with the binary. Planetesimals are immersed in 
a massive, axisymmetric gaseous disk, characterized by 
surface density E(r) specified in W2 2.11 

Assuming e « 1 the secular (averaged over the plan- 
etesimal and binary orbital motion) disturbing function 
for a planetesimal with semimajor axis a and eccentricity 
vector e = (k,h) = (ecosw, esinzu) (with apsidal angle 
w counted from the binary apsidal line, which is assumed 
fixed 2 ) is (Murray & Dermott 1999) 
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is the precession frequency of planetesimal orbit due to 
the disk potential U d . Here n b = [G(M P + M s )/ag] 1/2 
and n = (GMp/a 3 ) 1 ^ 2 are the mean motions of the bi- 
nary and planetesimal, respectively. The contribution to 
R proportional to vo d arises from expansion of the disk 
potential along the eccentric planetesimal orbit and av- 
eraging over its mean longitude. 
Evolution equations for h and k are written using 

dh/dt = (na 2 )' 1 dR/dk, dk/dt = - (na 2 )' 1 dR/dh as 



dh 
~db 



(A + zu d )k - B, 



dk 
~dt 



-(A + w d )h. (5) 



These equations agree with the work of Marzari & Scholl 
(2000) as long as w d = 0. 

We write down the solution for e(t) — ef rco (i) + 
eforced(i), where 
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Thus, free eccentricity vector ef ree rotates at a rate A + 
w d around the endpoint of the fixed forced eccentricity 
vector eforccd- Note that setting w d — we reproduce 
the solution obtained by Heppenheimer (1978). 

Planetesimals starting on circular orbits have ef rco — 
Cforced so that their eccentricity oscillates with amplitude 
2ef orce d and period 2ir/(A + vu d ). 

2 Precession period of the secondary orbit exceeds important 
timescales of the problem. 



2.1. Disk model. 

We model the disk as a constant M disk extending out 
to the outer truncation radius r t . Numerical simulations 
of accretion disks in binaries suggest that r t ~ (0.2 — 
0A)ab (Zhou et al. 2012), depending on e& and \i. In our 
study we will commonly take r t = 0.25a;,. 

Constant M assumption is a necessary simplification, 
which ignores the details of the disk structure at r ~ r t . 
Assuming viscosity v in the disk to be well described by 
the radius-independent effective a-parameter (Shakura 
& Sunyaev 1973), one can write 



M ClM 

E W = q = O 2 ' 

Ativ oiraci 



(8) 



From this it is obvious that if the disk temperature scales 
as Tex r~ q then the surface density behaves as 
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where tq is some fiducial radius and Eo = E(ro). 

Models of protoplanetary disks typically find q to be 
close to 1/2, so that p « 1. In particular, the passive 
disk model of Chiang & Goldreich (1997) has q — 3/7, so 
that p — 1 = 1/14. Outer parts of a disk in a binary are 
additionally heated by the radiation of the companion 
and tidal dissipation, so that q may be lower than 3/7 
even for a passive disk. For simplicity, in our calculations 
we will assume p = 1, which corresponds to a classical 
Mestel disk (Mestel 1963) if r t -)• oo. 

Assuming that the disk has power law profile © with 
p = 1 all the way to r t we can express its surface density 
via the total disk mass M d enclosed within r t as 
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where r t .5 = r t /(5AU), and ri = r/(lAU). Interestingly, 
gas surface density at 1 AU in such a disk with M d = 
O.OIMq is not very different from that in a Minimum 
Mass Solar Nebula (MMSN; Hayashi 1981). 

2.2. Precession due to the disk. 

A disk with the density profile ([9]) with p — 1 extending 
to infinity is known to have constant circular velocity 
(Mestel 1963) 
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Expressing dU d /dr from this relation and plugging it in 
equation Q we find vo d = —TtGT,(r)/(nr), where from 
now on we use r instead of a. Even though the cir- 
cumprimary disk in our problem is truncated at rt this 
expression should still be able to give us a reasonable es- 
timate of the precession rate w d due to the disk potential 
for r <r t . Using equation (|10[) we find 
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Note that za d varies rather weakly with r, as r 1 I 2 , which 
is consistent with Batygin et al. (2011). 
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2.3. Planetesimal eccentricities. 

To assess the role of disk-driven precession on secular 
evolution of planetesimals we evaluate 
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where a^o = a&/(20AU). In making this estimate we 
neglected the term quadratic in e b in equation @. It is 
obvious that dynamics of planetesimals at several AU is 
strongly affected by the disk-driven precession. Indeed, 
\w d \ exceeds A for 

4.6Auf^^Y /2 , (14) 
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i.e. over essentially the whole assumed extent of the disk 
even for the disk mass as small as <~ 10 -2 M Q . Thus, 
if we are interested in planet formation at 2 — 3 AU we 
can neglect planetesimal precession due to the secondary 
compared to the disk-driven precession, i.e. neglect A 
compared to vj d in equation ([7]) and other formulae. 

Equation ([7]) then predicts that the amplitude of ec- 
centricity oscillations is 

e - W = S»He*lf? ,15, 
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where we again neglected e 2 term in equation ([3]). This 
is to be compared with 
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which one finds neglecting disk-driven precession, i.e. 
dropping w d in equation ([T)l. It is obvious that neglect- 
ing disk-driven precession leads to an overestimate of the 
planetesimal eccentricity at ~AU separations by more 
than an order of magnitude. This has important conse- 
quences for planetesimal growth as we discuss further. 

2.4. Gas drag. 

Equation (fT5")) accounts for the presence of the disk 
only through the precession caused by its gravity. How- 
ever, for small planetesimals the effect of gas drag is also 
important. Assuming quadratic drag force in the form 
F « —wpg/(pd) (here d and p are the object's radius 
and bulk density, p g » £//i is the gas density and h is the 
disk scale height) we account for its effect on planetesi- 
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mal dynamics by adding terms —D{h, k} (h 2 + k/ 
the first and second equations ([5]), respectively (Marzari 
& Scholl 2000). Here 

D = nPJ 7 =n -/v ^ 

pd pd n 

For small planetesimal sizes (to be specified later by 
equation (j2T)]) ). in the gas drag-dominated regime, the 
drag force balances eccentricity excitation due to the sec- 
ondary, i.e. the B term in the first equation (JS|). This re- 
sults in the following estimate for the gas drag-mediated 
planetesimal eccentricity: 
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This expression agrees with Paardekooper et al. (2008) 
and predicts that e gas °c d 1 ! 2 . As a result, for small 
bodies one finds e gas < e dlsk . 

The transition between the drag-dominated behavior 
(fT9)l and the drag-free eccentricity scaling (TTBl occurs 
at the planetesimal size d gas where these two equations 
yield the same eccentricity: 

4nrS 2 15 r M s M p r t r 

dgas ~ B hp etorccd ~ 8n eb hM d pa 4 b [ UJ 
e b r/h 0.01 M pA r tj5 ri 
~ 0.5 30 M d /M s p 3 <20 ' 

where p 3 = p/(3 g cm~ 3 ), M p l = M p /M and we have 
used equation pU|) . 

Planetesimal eccentricity behaves according to formula 
(fT9]) for d < dgas, and switches to drag-free regime f(T5|) 
for d > d g as , see Figure [TJ The dependence of d gas on 
gas disk density and mass — d gas oc M^ — is some- 
what counter-intuitive, since higher gas density results 
in stronger drag, making it more important for larger 
bodies. However, the gas-free planetesimal eccentricity 
(fTS"|) is itself a function of M d and decreases faster with 
increasing M d than does e gas , explaining the nontrivial 
d g as(M d ) dependence. 

3. IMPLICATIONS FOR PLANET FORMATION. 

Planetesimals grow in mutual collisions as long as their 
encounter velocity v co u (measured at infinity) is such that 
collisions do not result in the net loss of mass. The 
conditions for this depend, in particular, on planetesi- 
mal size and on whether planetesimals are strength- or 
gravity-dominated. Using results of Leinhardt & Stewart 
(2012) for collisions of equal mass (the most disruptive) 
strengthless bodies we roughly estimate the condition for 
planetesimal growth to be 
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where the escape speed from the surface of an object of 
radius d and bulk density p is 

1/2 



8 TV 



Gp 



1 „ _i 1/2 , 

1.3 m s p 3 d\ 



(23) 



(here d\ = d/(lkm)). It becomes harder to break plan- 
etesimals when they are small enough for their internal 
strength to dominate over the gravitational energy, which 
is expected to happen for d < d s ~ 10 km (Holsapple 
1994). 

Planetesimal collisions occur at velocity of order 
Vcoiiif) ~ e dmk VK, where vk = nr is the Keplerian speed. 
Using expression (fTS")) we find 
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Plugging equations (|23|) and ([M} into the condition (|22|) 
we find that erosion in equal-mass planetesimal collisions 
is avoided for bodies with d > d co u , where 
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Thus, for the fiducial binary parameters adopted here 
and for M d ~ 10~ 2 M Q only planetesimals larger than 
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« 35 km would be able to grow at 1 AU. At 2 AU - 
the semi-major axis of 7 Cephei Ab — only bodies larger 
than 200 km in radius would be able to survive in equal- 
mass collisions. 

However, in the absence of a disk the problem is much 
worse: evaluating collisional velocity as v co u = e D ^ dlsk Vx 
using equation (fTT)) and applying condition (|2"2"|) one finds 
than in the absence of disk-induced precession only plan- 
etesimals larger than 
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are able to survive in equal-mass collisions. Clearly, col- 
lisional barrier appears far more severe if one disregards 
the effects of the disk on the secular evolution of plan- 
etesimals. We compare the behavior of d gas , d co u, and 



,n/disk 

oil 



as a function of r in Figure [TJ 



We also point out that d co ii is very sensitive to the 
binary semi-major axis ab, unlike d^^ lsk , see equations 
dM} and (|2"6]): increasing a b from 20 AU to 30 AU re- 
duces dcou by a factor of 5. The size of the region where 
disk-driven precession dominates secular evolution also 
expands rapidly with increasing a&, see equation ([T4|) . 
To summarize, properly accounting for the disk gravity 
considerably alleviates the collisional barrier in binaries, 
certainly for r < 1 AU. 

4. PLANETESIMAL ACCRETION IS POSSIBLE IN MASSIVE DISKS 

We now propose a solution to the problem of planetes- 
imal accumulation in binaries, raised in [JTJ We argue 
that if 




• disk is massive, Md ~ 0.1 Af, 



Fig. 1. — Characteristic planetesimal sizes vs. radius for two 
disk masses Mj: (a) O.OIMq and (b) O.IMq. We display d co u 

(solid, eq. |25|). d ga s (dot-dashed, eq. |20| ), d^ f^ (short- dashed, 
eq. 1261 ). and planetesimal radius d 3 = 10 km below which we con- 
sider objects as strength-dominated (long- dashed). The two latter 
sizes do not depend on M d . Calculations are done for e b = 0.4, 
M a = 0.4Mq , a b = 20 AU, M p = Mq (typical for small separation 
binaries, Chauvin et al. 2011), rt = 5 AU, and r/h = 30. Plan- 
etesimals in the shaded region ("danger zo ne" ) get destroyed in 
equal-mass collisions according to criterion J22H precluding plane- 
tary growth at corresponding separations. Accretion-friendly zone 
is to the left of the vertical dotted line in each plot; it is wider for 
higher M d and extends to Fa 2.5 AU for M d = O.IMq. 
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• planetesimals are strength- dominated below 
km, 
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then the fragmentation barrier can be overcome even at 
separations of w 2 AU, where planets in several binaries 
are found. 

Equation (|2"5"|) shows that higher Md results in smaller 
planetesimal size d co u, below which strengthless objects 
are destroyed or eroded in equal-mass collisions. Smaller 
planetesimals are (1) more resistive to collisional erosion 
because of their internal strength and (2) stronger af- 
fected by gas drag. When the latter dominates, plan- 
etesimal velocities are reduced and collisions are less de- 
structive. However, increasing Md reduces not only d co u 
but also dg as in such a way that d co u/d gas stays con- 
stant. At 1 AU this ratio is about 30 so that independent 
of Md there is still a significant "danger zone" between 
the planetesimal size d gas below which gas drag lowers 
v co u helping accretion and the radius d co u above which 
colliding strengthless bodies can grow, see Figure [1] 

On the other hand, equation (J2HJ) predicts that d co u < 
d s at r = 2 AU for M d /M s = 0.2, if internal strength 
dominates over the gravitational binding energy of the 
body with d s — 10 km making it harder to erode or 
destroy. This is a resolution of the collisional barrier 
problem in binaries that we favor in this work. 

To avoid fragmentation barrier we thus require that 
dcoii < d s . Using equation (|25[) we can rephrase this 
condition in the form of a lower limit on the disk mass 



at a given separation a p i from the primary: 
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In Figure [5] we illustrate this constraint as a function of 
the binary semi-major axis, for different values of a p i. 
It is clear that in very small separation binaries with 
ab = 10 AU growing planets even at 1 AU requires a 
massive disk, M p rs 0.2M s . 

External companions in giant planet-hosting binaries 
typically have mass M s « 0.4M q (Chauvin et al. 2011), 
meaning that our scenario of planet formation at 2 AU 
needs Md ~ O.IMq. Such disk mass may seem high but 
it is also the case that planet-hosting binary systems with 
ab ~ 20 AU that we consider here contain more mass in 
total than the descendants of the typical T Tauri stars. 

One might worry that such massive disks would be 
prone to gravitational instability (GI). With the density 
profile (fT0|) we estimate the Toomre Q = nc s /(nGYi) (c s 
is the sound speed) as 
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Thus, even for M d = 0.1M p w 0.1M Q the disk is at 
most marginally unstable to GI at 2 AU. However, even 
if it were unstable, the surface density and optical depth 
at this distance would be so high that the cooling time 
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Fig. 2. — Plot of ai, — Md/M s phase space illustrating conditions 
under which planets can form in binaries, assuming d s = 10 km. 
Different curves show the relation H27H for different values of the 
planetary semi-major axis a p i, interior to which planet formation 
is possible: a p i = 1 AU (solid), a p i = 2 AU (dotted), a p \ = 3 
AU (dashed). Calculation assumes = 0.4, fixed rt/aj, = 0.25, 
at, = 20 AU, and M p = Mq. At a given a p i fragmentation barrier 
is avoided and planet formation proceeds smoothly through the 
planetesimal stage to the right (and above) of the corresponding 
line. 



would far exceed the local dynamical time, making planet 
formation by direct disk fragmentation impossible (Gam- 
mic 2001; Rafikov 2005). Instead, the disk would slowly 
evolve under the action of gravitoturbulence (Rafikov 
2009). 

On the other hand, high M4 simplifies planet forma- 
tion in other ways. In particular, planets in these systems 
are quite massive (Chauvin et al. 2011) and larger Md 
provides mass reservoir for their assembly. Higher sur- 
face density of the protoplanetary disk also means larger 
isolation mass (Lissauer 1993) possibly making it high 
enough at 2 AU to trigger core accretion without the 
need to go through the long-lasting stage of giant impacts 
(Chambers 2004). Higher £ likely implies larger dead 
zone (Gammie 1996) in the disk, providing quiet condi- 
tions for planetesimal formation and growth, and result- 
ing in smaller viscosity, which, possibly, means longer 



disk lifetime. The timescale on which planets form also 
goes down as M<j increases. 

Another potential solution to the collisional barrier 
problem in binaries is the direct formation of large plan- 
etesimals by e.g. streaming and/or gravitational insta- 
bilities (Johansen et al. 2007; Thebault 2011). Large 
Md and small d co u are helpful for this mechanism as well 
since to overcome fragmentation barrier in a massive disk 
such instabilities would only need to produce bodies with 
sizes of tens of km, rather than ~ 10 3 km dwarf planets. 

5. DISCUSSION. 

We now mention several additional factors that may 
strengthen or weaken our conclusions. First, our colli- 
sional growth condition (|22|) may be too stringent. Pre- 
viously, using a more refined fragmentation criterion 
Thebault (2011) found that in HD 196885 (M p = 1.3M , 
ab = 21 AU, ef, = 0.42) planetesimal growth is possible 
even in the absence of disk-driven precession at 2.6 AU as 
long as the planetesimal size exceeds 250 km. However, 
according to our formula (j2"oT) with the same assumptions 
(same system parameters and Wd = 0), growth is possi- 
ble only for d > 10 3 km. Thus, our fragmentation cri- 
terion likely overestimates planetesimal size above which 
objects grow efficiently, and, in fact, it might be easier io 
overcome the fragmentation barrier with the more real- 
istic growth condition than our simple criterion (f22|) . At 
a given distance this would lower the value of Md needed 
to overcome fragmentation barrier. 

Second, even if planetesimals are collisionally weak 
their growth may still proceed mainly via unequal-mass 
collisions (which more frequently result in mergers) if the 
number of relatively massive objects is small (Thebault 
2011). Outward migration of planets by scattering of 
planetesimals has also been invoked (Payne et al. 2009) 
to explain planets on AU-scale orbits in small separation 
binaries. 

On the other hand, there are also factors complicat- 
ing planetesimal growth. In particular, eccentricity of 
the gaseous disk induced by the companion may affect 
planetary growth at small sizes (d < d gas )- Also, gas 
drag-induced inspiral of planetesimals may deplete the 
disk of some solids. The relative importance of these 
factors for planet formation in binaries will be assessed 
in the future. 
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